Gravitational spectrum of black holes in the Einstein- Aether theory. 
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Evolution of gravitational perturbations, both in time and frequency domains, is considered for a 
spherically symmetric black hole in the non-reduced Einstein- Aether theory. It is shown that real 
oscillation frequency and damping rate are larger for the Einstein-Aether black hole than for the 
Schwarzschild black hole. This may provide an opportunity to observe aether in the forthcoming 
experiments with new generation of gravitational antennas. 

PACS numbers: 04.30.Nk,04.50.+h 



o 
o 

(N 
C 

Hi 

in 



(N 

> 

(N 
(N 



o . 

^: 
^— > , 



X 



One of the most intriguing issues of modern physics 
consists in attempts to go beyond local Lorentz symme- 
try pj. In theory of gravity, breaking of local Lorentz 
invariance leads to a general relativity coupled to a dy- 
namical time-like vector field u a , called "aether". More 
exactly, u a breaks local boost invariance, while rotational 
symmetry in a preferred frame is preserved [2J . Thereby, 
aether is a kind of locally preferred state of rest at each 
point of space-time due-to some unknown physics. Re- 
cently observable consequences of Einstein- Aether theory 
attracted considerable interest Q. Gravitational conse- 
quences of Local Lorentz symmetry violation must show 
themselves in radiative processes around black holes. It 
is known that gravitational radiation damping of binary 
pulsars orbits reproduces the weak field general relativity 
at lowest post- Newtonian order Q. Yet, the significant 
difference between Einstein and Einstein- Aether theories 
should be seen in the regime of strong field, for instance 
in observing of the characteristic quasi-normal spectrum 
of black holes. Thus, existence of aether could be tested 
in the forthcoming experiments with new generation of 
gravitational antennas. Motivated by the above reasons, 
in a previous letter || we developed a method for finding 
of the quasinormal modes for the perturbations of met- 
rics which are not known analytically, but instead are 
given only numerically in some region near black holes. 
That is the case of the Einstein- Aether black holes found 
in Q . In Q , there were found the quasinormal modes for 
test scalar and electromagnetic fields in the vicinity of the 
Einstein- Aether black holes. It was shown that the scalar 
and electromagnetic quasinormal modes in the Einstein- 
Aether theory, have larger real oscillation frequency and 
damping rate than those of the Schwarzschild black holes 
in the Einstein theory. As quasinormal spectrum does 
not depend on the spin of the field in cikonal regime, 
qualitatively the same QN behavior was suggested in [f| 
for the gravitational perturbations as for scalar and elec- 
tromagnetic ones. In the present work we show that it 
is indeed true, and analyze gravitational perturbations 
of the Einstein- Aether black holes both in frequency and 
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time domain. 

We shall start from the lagrangian of the full Einstein- 
Aether theory forms the most general diffeomorphism in- 
variant action of the space-time metric g a b and the aether 
field u a involving no more than two derivatives given by 



L=-R-K 



ab 



y a u m V b u n - \{g ab u a u b - 1), (1) 



here R is the Ricci scalar, A is a Lagrange multiplier 
which provides the unit time- like constraint, 



K 



ab 



where the Ci are dimensionless constants. 

Spherically symmetry allows to fix C4 = 0. In this 
letter, following [||, we shall consider the so-called non- 
reduced Einstein- Aether theory, for which C3 = 0, and we 
can use the field redefinition that fixes the coefficient ci 



(>2 



2 - 4ci + 3ci 



so that ci is the free parameter. 

The metric for a spherically symmetric static black 
holes in Eddington-Finkelstein coordinates can be writ- 
ten in the form IfJ: 



ds 2 



N{r)dv 2 - 2B{r)dvdr - r 2 dn 2 , 



(2) 



where the functions N(r) and B{r) are given by numeri- 
cal integration near the black hole event horizon 0. One 
can re-write this metric in a Schwarzschild like form: 



ds 2 = -N{r)dt 2 + 



B 2 {r) 
N{r) 



dr 2 



2 dfl 2 



(3) 



Since the background value of aether coupling, de- 
termined by constants c$ is small in comparison with 
the background characteristics of large black hole, de- 
termined by the mass of the black hole M. There- 
fore, the background black hole metric is, in fact, the 
Schwarzschild metric slightly corrected by the aether. 
Thus, one can neglect small perturbations of aether, 
keeping only linear perturbations of Ricci tensor. Then 
the perturbation equations with unperturbed aether have 
the form: 



5R a p = 0. 



(4) 



2 



The general form of the perturbed metric, according to 
Chandrasckhar designations, is 



ds 2 



e 2v di 2 - e 2 ^{d$ 2 - u;dt - q 2 dr - q 3 ddf 



- 2fl2 dr 2 



Here 



N(r), 



- 2 ^d0 2 



N(r)/B 2 (r), 



r 2 sin 2 i 



(5) 



(6) 



Let us introduce new variables 

Qik = Qi,k - Qk,i-> Qio = 1i,o - i,k = 2,3. (7) 

Here we used 0,1,2,3 for t, (f>, r and 9 coordinates re- 
spectively. Here we shall consider the axial type of gravi- 
tational perturbations. The equations (11), (12) on p. 143 
of Q can be reduced to a single equation 



N dQ 
B \Br 2 dr 

where we used 



r 2 d 2 Q , 3 J 



B dt 2 



de 



1 dQ 

in 3 e oe 



= 
(8) 

(9) 



The following representation of the function Q(t, r, 9) 

Q(t, r, 9) = rRi{t, r)C^{ 2 {9) (10) 

leads to separation of angular variable 9. 

Finally we have the wave-like equation for the radial 
coordinate 



d 2 * 



+ [u? - F(r))# = 0, 



ilr , = - — \ dr. (11) 



dr 2 ' v ~ ' v "~ N(r) 
with the effective potential 

(£ + 2){£-l) 2N 2 (r) 1 d{N{r)/B{r)) 



V{r) = N(r)- 



B 2 (r)i 



dr* 



(12) 

Quasi-normal modes of asymptotically AdS black holes 
have been studies recent years extensively, because of 
their interpretation in Conformal Field Theory Q with 
some specific boundary conditions. In astrophysically 
relevant problem, one should require natural boundary 
conditions for QN modes of purely in-going waves at the 
event horizon and purely out-going waves at spatial in- 
finity 



±00. 



(13) 



Under these boundary conditions, the quasinormal 
modes were studied in a great number of papers [9(, yet 
in those cases the background metric and the effective 
potential were known in analytical form. For the case of 



Einstein- Aether theory, we are in position to apply the 
method developed in our previous paper Here we 
shall give only a brief summary of the whole procedure 
of@. 

We approximate the numerical data for the metric by 
a fit of the form 



N(r) 



N N 

E 

i=0 



B(r) 



N B 

E 

i=0 



A i. 



i+E 6 



i+E & r v 



i=i 



which are substituted into equations (fTTj) and (TT2"|) . The 
numbers Nn and Nb determine the number of terms in 
the polynomials and are chosen in order to provide best 
convergence of the WKB series. Coefficients a[ N \ b\ N \ 



AB) r(B) 



are determined by the fitting procedure. The 



WKB expansion has the form 



^-E^ = 



(14) 



where the correction terms of the i-th WKB order A, can 
be found in [l(| [ll[ and [12], Q = V — ui 2 and Qq means 
the i-th derivative of Q at its maximum. 

Alternatively, we shall use the above mentioned fits 
of the metric functions, and consequently of the effective 
potential, in the time-domain analysis: using the integra- 
tion scheme described for instance in [l3 ]. In detail, we 
used a numerical characteristic integration scheme, based 
in the light-cone variables u = t — r+ and v = t + r± . In 
the characteristic initial value problem, initial data are 
specified on the two null surfaces u = u Q and v — v . 
The discretization scheme applied, is 



*(JV) = $(W) + %(E) - 

*(W) + V(E) 



A 2 V(S)- 



+ 0(A 4 ) 



(15) 



where we have used the definitions for the points: N — 
{u + A, v + A), W = (u + A, v), E = (u, v + A) and 
S = (u, v). 

The application of the above two methods shows ex- 
cellent agreement: for instance the fundamental mode 
0.7686 — 0.1887i in time domain is very close to the WKB 
value 0.769470 — 0.187783i for ci = 0.4, as can be seen in 
Fig. 2. From the obtained numerical data in Table I and 
time domain pictures in Fig. 1-2, one can see that when 
increasing c\, both real oscillation frequency and damp- 
ing rate are increasing. Even for a small aether c\ ~ 0.1, 
the increase in Reu and Imuj is of about half percent, 
and could, in principle, be detected by new generation 
of gravitational antennas. For larger ci, the difference 
between, Schwarzschild and Einstein-Aether QNMs can 
be very significant and reach six-seven percents. From 
the Table I, it is evident that both real and imaginary 
parts of u> grows when increasing the multipole number 



3 



Table I: Axial gravitational perturbations for the non-reduced Einstein- Aether theory: The fundamental mode. 



Cl 


1 = 2 


£ = 3 


£ = 4 


e = 5 


0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.77 


0.751958 - 0.179578; 
0.757144 - 0.180786; 
0.762871 - 0.184443; 
0.769470 - 0.187783; 
0.777059 - 0.192176i 
0.786784 - 0.198087; 
0.799391 - 0.206961; 
0.811083 - 0.216302; 


1.206905 - 0.187387z 
1.215669 - 0.188850i 
1.225598 - 0.192601; 
1.236966 - 0. 196172i 
1.250474 - 0.200794; 
1.267215 - 0.207061i 
1.289844 - 0.216264; 
1.311449 - 0.225807; 


1.629395 - 0.190378i 
1.641492 - 0.191932i 
1.655131 - 0.195746i 
1.670792 - 0.199412; 
1.689421 - 0.204151; 
1.712545 - 0.210574; 
1.743878 - 0.220000; 
1.773893 - 0.229784; 


2.038526 - 0.191824; 
2.053816 - 0.193422; 
2.071021 - 0.197268; 
2.090803 - 0.20098H 
2.114345 - 0.205777; 
2.143595 - 0.212276; 
2.183271 - 0.221813; 
2.221342 - 0.231720; 



£. Here we considered only axial gravitational perturba- 
tions, which are iso-spectral with polar gravitational per- 
turbations for Schwarzschild black holes. For black holes 
in Einstein- Aether theory this iso-spectrality will be bro- 
ken, and QNMs for polar perturbations should slightly 
differ from axial, when considering the full perturbations 
of Einstein- Aether equations. The same breaking of iso- 
spectrality happens, for instance, when perturbing dila- 
ton black holes or black holes in higher than four di- 
mensional space-times In our approach, the per- 
turbations of aether were neglected in comparison with 



Figure 1: Evolution of axial gravitational perturbations in 
time domain ci = 0.4, £ = 3, non-reduced theory (red line) in 
comparison with the Schwarzschild case (blue line). 
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perturbations of the metric of a large astrophysical black 
hole. Therefore this difference between axial and polar 
QN spectra was neglected as well. 

Note that we used here the method based on the sup- 
position that QN frequencies are determined mainly near 
the peak of the potential barrier, while behavior of the 
potential barrier far from black hole is not significant. 
Even despite this idea was inspired by WKB approach, 
it is not dependent on WKB technique, as was shown 
here by computations in time domain. 

I 

Figure 2: Evolution of axial gravitational perturbations in 
time domain ci =0.1 (green line), ci = 0.4 (red line), I = 
2, non-reduced theory in comparison with the Schwarzschild 
case (blue line). 
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